Instabilities, nonhermiticity and exceptional points in the cranking model 
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A cranking harmonic oscillator model, widely used for the physics of fast rotating nuclei and 
Bose- Einstein condensates, is re-investigated in the context of PT-symmetry. The instability points 
of the model are identified as exceptional points. It is argued that - even though the Hamiltonian 
appears hermitian at first glance - it actually is not hermitian within the region of instability. 
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Quantum instabilities are attracting considerable at- 
tention in a variety of physical situations. They can be 
associated with the formation of solitons and vortices in 
Bose-Einstein condensates with a sudden change of 
the moment of inertia of a rotating nucleus (see, for exam- 
ple, ReflU and references therein) and a transition from 
one- to two-dimensional nuclear rotation [sj . A particular 
example of interest is the Hamiltonian 

H = huji{a\ai + -) + huj2{ala2 + -) + 

+ihgi{a\a2 — a^ai) — ihg2{a\al — a2ai) (1) 

used in condensed matter physics to describe in a sim- 
ple way the interaction between an atom and a radiative 
field |4j. Note that the bi- linear form of ^ corresponds 
to a linearised version of some more general interactions. 
As discussed below this may bring about an instability. 
Higher order terms may or may not remove such insta- 
bility. 

Using the standard relations [h — m — 1) 
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where Xi = x,X2 = y, and choosing special values for 
the strength constants gi = ^{uji + a;2)/2^wia'2 and 
52 = ^^{(^2 — wi)/2y/a}iW2, one recognises the well-known 
cranking Hamiltonian (Routhian) 
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which has been applied in nuclear physics [5L |6| and 
for rotating Bose-Einstein condensates (cf Reful). The 
Hamiltonian ([4]) appears as the sum of hermitian opera- 
tors and is expected - naively at first glance - to be itself 



a hermitian operator. The same holds when ([4]) is writ- 
ten in second quantised form ([T]). In the following we 
explore the formal character of the instability points of 
H and iJ. We argue that the operators are no longer 
hermitian at these points, in fact, we show, that these 
points are exceptional points (EP) [1, Q. 

Non-hermitian Hamilton operators have attracted 
widely spread interest during the recent years (see [lo|). 
be it in the context of effective theories 11 ill, or in the 
context of finding a hermitian equivalent [l2] or in the 
context of T'T-symmetry [l^ (PT is the product of the 
parity and time reversal operator). One specific aspect of 
non-hermitian operators are the EPs, being singularities 
of spectrum and eigenfunctions. As such, they are usu- 
ally of particular physical significance. They have been 
discussed in a great variety of physical applications: in 
optics 14 1, in mechanics llSll. as coalescing resonances 
1^ . 17 1, in atomic phvsics [iSj. and in more theoretical 
context in PT-symmetric models [l9| or in considering 
their mutual influence (20l | , to name just a few. In its sim- 
plest case they give rise to level repulsion being the more 
pronounced the nearer they lie to the real axis. Depend- 
ing on the particular situation they can signal a phase 
transition [21|. In the present case the EP is associated 
with the onset of an instability. Note that the Hamilto- 
nian ([4]) is symmetric under T'T-operation irrespective of 
a special choice of parameters (57 — > — J7 under T) ^2?|. 

It is well known that a Bogoliubov transformation of 
the Hamiltonian (j4]) 
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yields the form (cf [ 
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with the eigenmode energies 

(7) 

It is also known [2^ that becomes negative when 
the rotational speed lies between m\n[ujx,i^y) and 
uYa,yi{LOx,ijJy). In the following we assume that uj^ > LOy. 
At the points where the two eigenmodes vanish, that is 
when LO-^i = +w_ and a;_^2 = — W- coalesce, the matrix 
S in ([5]) becomes singular. This happens at the critical 
points rici = or f2c2 = signalling an instability. 

The coalescence is reminiscent of the behaviour of an 
EP. To confirm that we are in fact encountering a genuine 
EP and not a usual degeneracy, we have to analyse the 
eigenfunctions of the respective Hamiltonians. Of course, 
this is closely related to the singular behaviour of the 
Bogoliubov transformation B. 

To illuminate both, the underlying physics and the 
mathematical structure, it is convenient to construct the 
matrix U connecting the original canonical coordinates p 
and f with the quasi-boson operators qt and q\, that is 



(8) 



As a first step we aim at the generalised classical nor- 
mal mode coordinates P = (P+, P_) and X = (^+, X^) 
in which H assumes the form 
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This is achieved by solving the classical equations of mo- 
tion 
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which can be written in matrix form 
d 
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Note that (|4]) can be written as 

r) H 



H = 



(15) 



The solution of (|12p is obtained by exponentiation and 
reads 



(16) 



where T) — diag(— — ici;_, ia;_, icj-j.) is the diagonal 
form of = UW containing the eigenmodes. The 
columns of U and V are the right hand and left hand 
eigenvectors, respectively, of M.. Note that the eigen- 
modes are obtained from the non-symmetric matrix A^; 
from this classical view point it is therefore no surprise 
that some of the eigenvalues occurring in ^ may be 
complex. As the column vectors of lA and V form a bi- 
orthogonal system, we can choose V — U^^ . Also, we 
observe from the special form of p3)) that 
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While the explicit form of P(t) and X(i) is of little 
interest the essential point here is the classical insta- 
bility occurring for negative values of uPl^ that is for 
LOy < ^ < LOx- In fact, the harmonic oscillator po- 
tential has the 'wrong' sign in ^ for the coordinates 
P-{t) and X^{t). From we read off the classical 
'run away' solution in this parameter range yielding the 
~ exp(|ci;_|t) behaviour for position and momentum. 
The corresponding quantum mechanical behaviour is dis- 
cussed below. 

Using the form (|15p we aim at a form corresponding 
to ([6]), viz. 
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From ([5]) this implies that U must be normalised such 
that 7i = U^HqmU. The explicit form of the matrix el- 
ements of U are given in :23| , however the the quoted pa- 
per focuses upon significantly smaller values of f2 than the 
range of instability. The analytic form allows pertinent 
statements in general, and in particular an expansion in 
around the critical points fid = ujy and r2c2 = ^x- 
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FIG. 1: Relevant spectrum ±iO- as a function of in arbi- 
trary units. Parameters chosen are ujx = 3 and uiy — 2. The 
dashed lines indicate the imaginary part. 



The essential results are as follows: 

(i) When ^ [u!y,u!x] the commutators [^j,?]^] = Sj,k 
follow from [r,„,p„] = i5m,n- It guarantees that the bo- 
son operators are creation and annihilation operators for 
the excitations, in the present case, with energies w_ and 
ujj^. However, this holds only when uPi is positive; if cj^ 
is negative (fi € [ujy, uJx]) the commutator [q_, q^_] is neg- 
ative with the implication that the operators g_ and gl 
are no longer proper boson operators. We further note 
that V - as given by ([T7)) - no longer is the inverse of U 
for this parameter range. 

(ii) The end points of the instability region, i.e. the 
points fid = LOy and = can be clearly identified as 
EPs. In fact, while the two eigenvectors associated with 
the two distinct eigenvalues +uj- and — u;_ are obviously 
linearly independent, they become aligned, i.e. linearly 
dependent, at fici and figa where lj^ vanishes; this is the 
clear signature of an EP [24] . We recall: a genuine degen- 
eracy would have two linearly independent eigenvectors. 
EPs are a universal phenomenon occurring in spectra and 
eigenfunctions under variation of parameters. For hermi- 
tian operators they can occur only when such parameters 
are continued into the complex plane thus rendering the 
original hermitian operator effectively nonhermitian. 

EPs are square root singularities of the spectrum: in 
the present case the spectrum has a branch cut in fi rang- 
ing from ujy to ojx- When the eigenvalues +u}- and — 
are continued beyond the EP, they become imaginary (as 
was also noticed in f2^) for fi G [wyjWa;], again with op- 
posite sign (see Fig.l); clearly this contradicts H being 
hermitian for this parameter range. 

(iii) The correct normalisation enforced by ([8]) (to 
guarantee the correct commutation relations when ^ 
[LOy.LUx]) has the consequence that the leading terms of 
the components of the critical eigenvectors behave as 
(O — J7c)~^^^ when approaching the critical point. This 
particular singular behaviour - the forth root and the in- 
finity - is again a consequence of the eigenfunctions at an 
EP In fact, it has been shown in general [i^ that the 
scalar product of the two eigenfunctions - associated with 
the two coalescing levels - must vanish as a square root, in 
the present case as (fi — fic)^^^- As a consequence, when 



normalisation is enforced by dividing by the square root 
of the scalar product, the singular behaviour follows as 
indicated. Moreover, the forth root has the consequence 
that - for the wave function - a clockwise encircling of the 
EP in the fi-plane yields a result that has a phase that 
is different from that of a counterclockwise encirclement. 
In fact, considering yfz (taking z = 0. — fie), one obtains 
when z has described a full counterclockwise circle 
around zero and —i when going in the opposite direc- 
tion. This particular Riemann sheet structure has been 
experimentally established in microwave cavities 27|. It 
would be a challenge to confirm it in the present context 
with a BEC or with Raman scattering using an incident 
laser beam upon vibrational modes of a medium. 

So far, we have established the seemingly surpris- 
ing result that the Hamilton operator (|3|) - or its sec- 
ond quantised counterpart - fails to be hermitian when 
LOy < < bJx- The endpoints of this interval are EPs. 
We stress that this result is based on an analytic contin- 
uation obtained from the range fi < w^, or equally, from 
the range fi > w^,. These two (hermitian) ranges are of 
course also analytically connected. 

It appears apposite to contrast our findings with com- 
mon wisdom about the solutions of the Schroedinger 
equation of ([5]). In fact, if © is considered in isola- 
tion, the Hamiltonian appears perfectly hermitian, also 



for 



< 0. It has a continuous spectrum associated 



with the unbounded classical motion in the coordinates 
P- and The quantum mechanical wave function is 
asymptotically of the form exp(j|cj_ |X^/2) apart from a 
hypergeometric function. The crucial aspect explaining 
this apparent discrepancy lies in the transformation that 
brings us from ([3]) to ([5]). As long as fi ^ [wy,^^;] the 
two operators are equivalent up to a similarity transfor- 
mation. For fi G [wyjWa;] they are not. And the transfor- 
mation breaks down exactly at the EPs, the singularity 
that signals the instability point. 

This is a beautiful demonstration of a 'PT-symmetric 
operator yet with a special twist: (|4]) appears hermi- 
tian to the naked eye, but its spectrum is not real when 
fi € [a;y,Wa;]. While the operator is T^T-symmetric, the 
symmetry is broken by the state vector. Thus, in this pa- 
rameter range the hermitian form ^ is not its hermitian 
equivalent. 

In conclusion, we mention that the two exceptional 
points collapse into a diabolic point [28| when Uy = ujx] 
in this case fi = is a regular point with a genuine 
degeneracy for a)_ = 0. 
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